Project 2: Breakout Strategy¶
Instructions¶
Each problem consists of a function to implement and instructions on how to implement the function. The parts of the function that need to be implemented are marked with a # TODO comment. After implementing the function, run the cell to test it against the unit tests we've provided. For each problem, we provide one or more unit tests from our project_tests package. These unit tests won't tell you if your answer is correct, but will warn you of any major errors. Your code will be checked for the correct solution when you submit it to Udacity.
Packages¶
When you implement the functions, you'll only need to you use the packages you've used in the classroom, like Pandas and Numpy. These packages will be imported for you. We recommend you don't add any import statements, otherwise the grader might not be able to run your code.
The other packages that we're importing are helper, project_helper, and project_tests. These are custom packages built to help you solve the problems. The helper and project_helper module contains utility functions and graph functions. The project_tests contains the unit tests for all the problems.
Install Packages¶
In [1]:
import sys
!{sys.executable} -m pip install -r requirements.txt
Requirement already satisfied: colour==0.1.5 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 1)) (0.1.5)
Collecting cvxpy==1.0.3 (from -r requirements.txt (line 2))
Downloading https://files.pythonhosted.org/packages/a1/59/2613468ffbbe3a818934d06b81b9f4877fe054afbf4f99d2f43f398a0b34/cvxpy-1.0.3.tar.gz (880kB)
100% |████████████████████████████████| 880kB 17.4MB/s ta 0:00:01
Requirement already satisfied: cycler==0.10.0 in /opt/conda/lib/python3.6/site-packages/cycler-0.10.0-py3.6.egg (from -r requirements.txt (line 3)) (0.10.0)
Collecting numpy==1.13.3 (from -r requirements.txt (line 4))
Downloading https://files.pythonhosted.org/packages/57/a7/e3e6bd9d595125e1abbe162e323fd2d06f6f6683185294b79cd2cdb190d5/numpy-1.13.3-cp36-cp36m-manylinux1_x86_64.whl (17.0MB)
100% |████████████████████████████████| 17.0MB 2.7MB/s eta 0:00:01
Collecting pandas==0.21.1 (from -r requirements.txt (line 5))
Downloading https://files.pythonhosted.org/packages/3a/e1/6c514df670b887c77838ab856f57783c07e8760f2e3d5939203a39735e0e/pandas-0.21.1-cp36-cp36m-manylinux1_x86_64.whl (26.2MB)
100% |████████████████████████████████| 26.2MB 2.9MB/s eta 0:00:01
Collecting plotly==2.2.3 (from -r requirements.txt (line 6))
Downloading https://files.pythonhosted.org/packages/99/a6/8214b6564bf4ace9bec8a26e7f89832792be582c042c47c912d3201328a0/plotly-2.2.3.tar.gz (1.1MB)
100% |████████████████████████████████| 1.1MB 31.0MB/s ta 0:00:01
Requirement already satisfied: pyparsing==2.2.0 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 7)) (2.2.0)
Requirement already satisfied: python-dateutil==2.6.1 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 8)) (2.6.1)
Requirement already satisfied: pytz==2017.3 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 9)) (2017.3)
Requirement already satisfied: requests==2.18.4 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 10)) (2.18.4)
Collecting scipy==1.0.0 (from -r requirements.txt (line 11))
Downloading https://files.pythonhosted.org/packages/d8/5e/caa01ba7be11600b6a9d39265440d7b3be3d69206da887c42bef049521f2/scipy-1.0.0-cp36-cp36m-manylinux1_x86_64.whl (50.0MB)
100% |████████████████████████████████| 50.0MB 771kB/s eta 0:00:01
Requirement already satisfied: scikit-learn==0.19.1 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 12)) (0.19.1)
Requirement already satisfied: six==1.11.0 in /opt/conda/lib/python3.6/site-packages (from -r requirements.txt (line 13)) (1.11.0)
Collecting tqdm==4.19.5 (from -r requirements.txt (line 14))
Downloading https://files.pythonhosted.org/packages/71/3c/341b4fa23cb3abc335207dba057c790f3bb329f6757e1fcd5d347bcf8308/tqdm-4.19.5-py2.py3-none-any.whl (51kB)
100% |████████████████████████████████| 61kB 26.4MB/s ta 0:00:01
Collecting zipline==1.2.0 (from -r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/15/d3/689f2a940478b82ac57c751a40460598221fd82b0449a7a8f7eef47a3bcc/zipline-1.2.0.tar.gz (659kB)
100% |████████████████████████████████| 665kB 31.7MB/s ta 0:00:01
Collecting osqp (from cvxpy==1.0.3->-r requirements.txt (line 2))
Downloading https://files.pythonhosted.org/packages/ea/d4/50bf412bc9f3c4e84a7cfed014fe9c6d9b0b04b353a57325704905d1c176/osqp-0.6.3-cp36-cp36m-manylinux_2_5_x86_64.manylinux1_x86_64.manylinux_2_17_x86_64.manylinux2014_x86_64.whl (233kB)
100% |████████████████████████████████| 235kB 36.1MB/s ta 0:00:01
Collecting ecos>=2 (from cvxpy==1.0.3->-r requirements.txt (line 2))
Downloading https://files.pythonhosted.org/packages/04/da/aefd27c06a9179a7e5614d0d97c0384072d2d22800790690c661eb6f2f4a/ecos-2.0.12.tar.gz (142kB)
100% |████████████████████████████████| 143kB 37.5MB/s ta 0:00:01
Installing build dependencies ... done
Collecting scs>=1.1.3 (from cvxpy==1.0.3->-r requirements.txt (line 2))
Downloading https://files.pythonhosted.org/packages/b5/21/0391cc2f02ee80bbb44fec4f1f22d00066b690dc10148582f8421ce1f7bc/scs-3.2.3.tar.gz (1.2MB)
100% |████████████████████████████████| 1.2MB 23.0MB/s ta 0:00:01
Installing build dependencies ... done
Collecting multiprocess (from cvxpy==1.0.3->-r requirements.txt (line 2))
Downloading https://files.pythonhosted.org/packages/df/d0/c3011a5bb77f307e68682a5046cce1a2c6591267bf24b5bf3fc4130bb39d/multiprocess-0.70.12.2-py36-none-any.whl (106kB)
100% |████████████████████████████████| 112kB 36.4MB/s ta 0:00:01
Requirement already satisfied: fastcache in /opt/conda/lib/python3.6/site-packages (from cvxpy==1.0.3->-r requirements.txt (line 2)) (1.0.2)
Requirement already satisfied: toolz in /opt/conda/lib/python3.6/site-packages (from cvxpy==1.0.3->-r requirements.txt (line 2)) (0.8.2)
Requirement already satisfied: decorator>=4.0.6 in /opt/conda/lib/python3.6/site-packages (from plotly==2.2.3->-r requirements.txt (line 6)) (4.0.11)
Requirement already satisfied: nbformat>=4.2 in /opt/conda/lib/python3.6/site-packages (from plotly==2.2.3->-r requirements.txt (line 6)) (4.4.0)
Requirement already satisfied: chardet<3.1.0,>=3.0.2 in /opt/conda/lib/python3.6/site-packages (from requests==2.18.4->-r requirements.txt (line 10)) (3.0.4)
Requirement already satisfied: idna<2.7,>=2.5 in /opt/conda/lib/python3.6/site-packages (from requests==2.18.4->-r requirements.txt (line 10)) (2.6)
Requirement already satisfied: urllib3<1.23,>=1.21.1 in /opt/conda/lib/python3.6/site-packages (from requests==2.18.4->-r requirements.txt (line 10)) (1.22)
Requirement already satisfied: certifi>=2017.4.17 in /opt/conda/lib/python3.6/site-packages (from requests==2.18.4->-r requirements.txt (line 10)) (2023.5.7)
Requirement already satisfied: pip>=7.1.0 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (18.1)
Requirement already satisfied: setuptools>18.0 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (38.4.0)
Collecting Logbook>=0.12.5 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/2f/d9/16ac346f7c0102835814cc9e5b684aaadea101560bb932a2403bd26b2320/Logbook-1.5.3.tar.gz (85kB)
100% |████████████████████████████████| 92kB 31.8MB/s ta 0:00:01
Collecting requests-file>=1.4.1 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/77/86/cdb5e8eaed90796aa83a6d9f75cfbd37af553c47a291cd47bc410ef9bdb2/requests_file-1.5.1-py2.py3-none-any.whl
Collecting pandas-datareader<0.6,>=0.2.1 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/40/c5/cc720f531bbde0efeab940de400d0fcc95e87770a3abcd7f90d6d52a3302/pandas_datareader-0.5.0-py2.py3-none-any.whl (74kB)
100% |████████████████████████████████| 81kB 17.9MB/s ta 0:00:01
Requirement already satisfied: patsy>=0.4.0 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (0.4.1)
Requirement already satisfied: statsmodels>=0.6.1 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (0.8.0)
Requirement already satisfied: Cython>=0.25.2 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (0.29.7)
Collecting cyordereddict>=0.2.2 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/d1/1a/364cbfd927be1b743c7f0a985a7f1f7e8a51469619f9fefe4ee9240ba210/cyordereddict-1.0.0.tar.gz (138kB)
100% |████████████████████████████████| 143kB 38.5MB/s ta 0:00:01
Collecting bottleneck>=1.0.0 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/4e/b1/abd8986be57a9ed6f93ad685c55ebe3fe93ed924109a1978c1b659074a18/Bottleneck-1.3.7-cp36-cp36m-manylinux_2_5_x86_64.manylinux1_x86_64.manylinux_2_17_x86_64.manylinux2014_x86_64.whl (351kB)
100% |████████████████████████████████| 358kB 32.3MB/s ta 0:00:01
Collecting contextlib2>=0.4.0 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/76/56/6d6872f79d14c0cb02f1646cbb4592eef935857c0951a105874b7b62a0c3/contextlib2-21.6.0-py2.py3-none-any.whl
Requirement already satisfied: networkx<2.0,>=1.9.1 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (1.11)
Requirement already satisfied: numexpr>=2.6.1 in /opt/conda/lib/python3.6/site-packages (from zipline==1.2.0->-r requirements.txt (line 15)) (2.6.4)
Collecting bcolz<1,>=0.12.1 (from zipline==1.2.0->-r requirements.txt (line 15))
Downloading https://files.pythonhosted.org/packages/6c/8b/1ffa01f872cac36173c5eb95b58c01040d8d25f1b242c48577f4104cd3ab/bcolz-0.12.1.tar.gz (622kB)
100% |████████████████████████████████| 624kB 31.7MB/s ta 0:00:01
Complete output from command python setup.py egg_info:
Traceback (most recent call last):
File "<string>", line 1, in <module>
File "/tmp/pip-install-kitp0gzp/bcolz/setup.py", line 211, in <module>
cmdclass=LazyCommandClass(),
File "/opt/conda/lib/python3.6/site-packages/setuptools/__init__.py", line 129, in setup
return distutils.core.setup(**attrs)
File "/opt/conda/lib/python3.6/distutils/core.py", line 108, in setup
_setup_distribution = dist = klass(attrs)
File "/opt/conda/lib/python3.6/site-packages/setuptools/dist.py", line 333, in __init__
_Distribution.__init__(self, attrs)
File "/opt/conda/lib/python3.6/distutils/dist.py", line 281, in __init__
self.finalize_options()
File "/opt/conda/lib/python3.6/site-packages/setuptools/dist.py", line 476, in finalize_options
ep.load()(self, ep.name, value)
File "/opt/conda/lib/python3.6/site-packages/pkg_resources/__init__.py", line 2408, in load
return self.resolve()
File "/opt/conda/lib/python3.6/site-packages/pkg_resources/__init__.py", line 2414, in resolve
module = __import__(self.module_name, fromlist=['__name__'], level=0)
File "/tmp/pip-install-kitp0gzp/bcolz/.eggs/setuptools_scm-7.1.0-py3.6.egg/setuptools_scm/__init__.py", line 5
from __future__ import annotations
^
SyntaxError: future feature annotations is not defined
----------------------------------------
Command "python setup.py egg_info" failed with error code 1 in /tmp/pip-install-kitp0gzp/bcolz/
Load Packages¶
In [2]:
import pandas as pd
import numpy as np
import helper
import project_helper
import project_tests
Market Data¶
Load Data¶
While using real data will give you hands on experience, it's doesn't cover all the topics we try to condense in one project. We'll solve this by creating new stocks. We've create a scenario where companies mining Terbium are making huge profits. All the companies in this sector of the market are made up. They represent a sector with large growth that will be used for demonstration latter in this project.
In [3]:
df_original = pd.read_csv('../../data/project_2/eod-quotemedia.csv', parse_dates=['date'], index_col=False)
# Add TB sector to the market
df = df_original
df = pd.concat([df] + project_helper.generate_tb_sector(df[df['ticker'] == 'AAPL']['date']), ignore_index=True)
close = df.reset_index().pivot(index='date', columns='ticker', values='adj_close')
high = df.reset_index().pivot(index='date', columns='ticker', values='adj_high')
low = df.reset_index().pivot(index='date', columns='ticker', values='adj_low')
print('Loaded Data')
/opt/conda/lib/python3.6/site-packages/ipykernel_launcher.py:5: FutureWarning: Sorting because non-concatenation axis is not aligned. A future version of pandas will change to not sort by default. To accept the future behavior, pass 'sort=False'. To retain the current behavior and silence the warning, pass 'sort=True'.
Loaded Data
View Data¶
To see what one of these 2-d matrices looks like, let's take a look at the closing prices matrix.
In [4]:
close
Out[4]:
| ticker | A | AAL | AAP | AAPL | ABBV | ABC | ABT | ACN | ADBE | ADI | ... | XL | XLNX | XOM | XRAY | XRX | XYL | YUM | ZBH | ZION | ZTS |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| date | |||||||||||||||||||||
| 2013-07-01 | 29.99418563 | 16.17609308 | 81.13821681 | 53.10917319 | 34.92447839 | 50.86319750 | 31.42538772 | 64.69409505 | 46.23500000 | 39.91336014 | ... | 27.66879066 | 35.28892781 | 76.32080247 | 40.02387348 | 22.10666494 | 25.75338607 | 45.48038323 | 71.89882693 | 27.85858718 | 29.44789315 |
| 2013-07-02 | 29.65013670 | 15.81983388 | 80.72207258 | 54.31224742 | 35.42807578 | 50.69676639 | 31.27288084 | 64.71204071 | 46.03000000 | 39.86057632 | ... | 27.54228410 | 35.05903252 | 76.60816761 | 39.96552964 | 22.08273998 | 25.61367511 | 45.40266113 | 72.93417195 | 28.03893238 | 28.57244125 |
| 2013-07-03 | 29.70518453 | 16.12794994 | 81.23729877 | 54.61204262 | 35.44486235 | 50.93716689 | 30.72565028 | 65.21451912 | 46.42000000 | 40.18607651 | ... | 27.33445191 | 35.28008569 | 76.65042719 | 40.00442554 | 22.20236479 | 25.73475794 | 46.06329899 | 72.30145844 | 28.18131017 | 28.16838652 |
| 2013-07-05 | 30.43456826 | 16.21460758 | 81.82188233 | 54.17338125 | 35.85613355 | 51.37173702 | 31.32670680 | 66.07591068 | 47.00000000 | 40.65233352 | ... | 27.69589920 | 35.80177117 | 77.39419581 | 40.67537968 | 22.58516418 | 26.06075017 | 46.41304845 | 73.16424628 | 29.39626730 | 29.02459772 |
| 2013-07-08 | 30.52402098 | 16.31089385 | 82.95141667 | 53.86579916 | 36.66188936 | 52.03746147 | 31.76628544 | 66.82065546 | 46.62500000 | 40.25645492 | ... | 27.98505704 | 35.20050655 | 77.96892611 | 40.64620776 | 22.48946433 | 26.22840332 | 46.95062632 | 73.89282298 | 29.57661249 | 29.76536472 |
| 2013-07-09 | 30.68916447 | 16.71529618 | 82.43619048 | 54.81320389 | 36.35973093 | 51.69535307 | 31.16522893 | 66.48866080 | 47.26000000 | 40.69632003 | ... | 28.31939579 | 35.50113886 | 78.89018496 | 40.80179133 | 22.48946433 | 26.58233774 | 47.28094525 | 73.70108798 | 28.91218282 | 29.80384612 |
| 2013-07-10 | 31.17771395 | 16.53235227 | 81.99032166 | 54.60295791 | 36.85493502 | 52.28710814 | 31.16522893 | 66.71298151 | 47.25000000 | 41.10979324 | ... | 27.95794850 | 36.39419366 | 78.45068533 | 40.71427558 | 22.96796358 | 26.98284247 | 47.08340158 | 74.00785631 | 28.32368796 | 29.86156823 |
| 2013-07-11 | 31.45983407 | 16.72492481 | 82.00022986 | 55.45406479 | 37.08155384 | 53.72026495 | 31.85599537 | 67.47567196 | 47.99000000 | 42.22705062 | ... | 28.50011944 | 37.00430040 | 78.83102155 | 41.01571874 | 23.23113816 | 27.03872686 | 46.54333492 | 74.93774876 | 27.84909533 | 29.74612402 |
| 2013-07-12 | 31.48047700 | 16.90786872 | 81.91105609 | 55.35309481 | 38.15724076 | 53.98840397 | 31.81096287 | 67.76280247 | 48.39000000 | 42.53495620 | ... | 28.92482002 | 38.00346072 | 78.94089646 | 40.83096325 | 23.49431274 | 27.08529718 | 45.96422730 | 75.68549560 | 28.44708204 | 30.15979909 |
| 2013-07-15 | 31.72819223 | 17.10044125 | 82.61453801 | 55.47379158 | 37.79303181 | 53.84971137 | 31.95506689 | 68.41781897 | 48.12000000 | 42.57894271 | ... | 29.27723113 | 38.17146113 | 78.81411772 | 40.84068723 | 23.54216266 | 27.06666905 | 46.69299195 | 76.27027369 | 28.77929688 | 30.38106716 |
| 2013-07-16 | 31.59057266 | 17.28338516 | 81.62371841 | 55.83133953 | 37.10696377 | 53.88669607 | 32.15320992 | 67.55642741 | 47.48500000 | 42.68451033 | ... | 29.04229039 | 38.27314559 | 78.85637730 | 40.86013517 | 23.27898808 | 26.61959399 | 46.56936223 | 76.81670381 | 28.06740794 | 29.97701243 |
| 2013-07-17 | 31.38414330 | 17.76481650 | 80.74188897 | 55.84626440 | 37.23401341 | 54.06237335 | 32.26128793 | 67.43978064 | 48.04000000 | 42.80767257 | ... | 29.18686931 | 38.48977769 | 78.99160796 | 40.93792696 | 23.18328823 | 26.66616431 | 46.45874617 | 78.30261578 | 28.06740794 | 29.81346647 |
| 2013-07-18 | 31.58369168 | 17.73593062 | 81.74261676 | 56.03418797 | 37.53893253 | 53.91443458 | 32.15320992 | 67.69101984 | 48.19000000 | 42.52615889 | ... | 29.55735279 | 40.52346684 | 79.76918424 | 41.22964615 | 23.49431274 | 26.94558622 | 46.97929234 | 78.81069986 | 28.77929688 | 29.64992051 |
| 2013-07-19 | 31.79012104 | 17.55298671 | 81.45527908 | 55.15063572 | 37.70833205 | 54.37674323 | 32.30632044 | 67.49361761 | 48.07000000 | 42.20945601 | ... | 29.71096789 | 40.54999322 | 80.43688561 | 41.24909410 | 23.20721320 | 26.81518933 | 46.90121042 | 81.16898043 | 28.99760949 | 29.09194018 |
| 2013-07-22 | 32.20297975 | 17.47595770 | 81.99032166 | 55.32713852 | 38.08948096 | 54.54317435 | 32.24327493 | 67.29621538 | 48.28000000 | 42.17426681 | ... | 29.84651063 | 40.59420386 | 80.14952046 | 41.49219343 | 23.47038778 | 26.88970184 | 46.50429396 | 81.02518181 | 29.27287321 | 29.12080123 |
| 2013-07-23 | 31.97590746 | 17.37967143 | 81.94078068 | 54.37713815 | 37.53046256 | 53.28569482 | 33.03584705 | 66.62325323 | 48.07000000 | 42.56134810 | ... | 29.13265221 | 40.52346684 | 80.46224136 | 41.32688588 | 23.42253785 | 26.74067682 | 45.82758393 | 81.00601167 | 28.38063907 | 28.91877387 |
| 2013-07-24 | 32.17545584 | 17.81295964 | 80.78152175 | 57.17003539 | 36.96297418 | 52.49052395 | 32.82869752 | 66.14769330 | 47.80000000 | 42.42938857 | ... | 28.73506019 | 40.24051879 | 80.28475112 | 41.15185437 | 23.51823770 | 26.62890805 | 46.49128030 | 80.56503316 | 28.53250871 | 28.76484826 |
| 2013-07-25 | 32.10664605 | 18.13070432 | 81.46518728 | 56.90917464 | 37.47117273 | 53.26720248 | 32.94578204 | 65.62726924 | 47.79000000 | 42.88684829 | ... | 29.02421802 | 41.05399445 | 80.26784729 | 40.91847901 | 23.44646282 | 26.85244558 | 46.91422407 | 79.47217195 | 28.19080202 | 29.36130999 |
| 2013-07-26 | 31.37726233 | 18.38104862 | 81.88133151 | 57.23233050 | 37.93702140 | 54.06237335 | 33.12591207 | 65.60932358 | 47.64000000 | 42.71969954 | ... | 29.11457985 | 40.83294128 | 80.11571280 | 40.98654682 | 23.18328823 | 26.70342056 | 48.15052124 | 80.98684153 | 28.04842424 | 29.27472684 |
| 2013-07-29 | 31.19835688 | 18.51584940 | 81.57417743 | 58.11484449 | 38.16571074 | 53.98840397 | 33.08988606 | 64.93636143 | 47.17000000 | 42.66691573 | ... | 28.92482002 | 40.38199282 | 79.47336717 | 40.93792696 | 23.08758839 | 26.50782523 | 47.83819353 | 80.16240164 | 27.71620940 | 28.94763492 |
| 2013-07-30 | 30.86118893 | 18.48696352 | 81.43546269 | 58.83253602 | 37.86079161 | 53.84046520 | 33.21597708 | 66.16563896 | 47.36000000 | 43.04519972 | ... | 28.38264907 | 40.88599404 | 79.28742502 | 41.08378656 | 23.06366342 | 23.78811863 | 47.53237266 | 79.55844670 | 27.77316051 | 28.96206545 |
| 2013-07-31 | 30.77861719 | 18.63139292 | 81.73270857 | 58.73000866 | 38.52144972 | 53.87744989 | 32.99081455 | 66.22844876 | 47.28000000 | 43.44107832 | ... | 28.32843198 | 41.28388974 | 79.23671352 | 41.69639686 | 23.20721320 | 23.21996075 | 47.44778390 | 80.02819050 | 28.13385091 | 28.74031861 |
| 2013-08-01 | 31.68002538 | 18.66027880 | 82.66407899 | 59.26808264 | 38.32664028 | 54.36749706 | 33.17995107 | 67.16162295 | 47.70000000 | 43.93372725 | ... | 28.97000093 | 41.69062757 | 78.37461808 | 41.71584481 | 23.70963740 | 23.46212640 | 48.08545297 | 80.97725310 | 28.61793539 | 29.07775945 |
| 2013-08-02 | 31.91397865 | 18.21736196 | 82.70371177 | 60.02912118 | 38.38593011 | 54.07161953 | 33.09889256 | 66.92832940 | 47.45000000 | 43.87214613 | ... | 28.87060292 | 41.12473146 | 77.71536862 | 41.78391262 | 23.92496206 | 23.53663891 | 48.40428750 | 80.52668616 | 28.61793539 | 29.82977047 |
| 2013-08-05 | 31.61121560 | 18.45807764 | 82.64426260 | 60.92591114 | 37.86926159 | 54.58015904 | 32.82869752 | 66.60530757 | 47.63000000 | 43.61702437 | ... | 28.60855363 | 41.00978381 | 77.41109964 | 41.52136535 | 24.09243679 | 23.54595298 | 48.68408107 | 80.46916901 | 28.39962278 | 30.12864664 |
| 2013-08-06 | 31.70754930 | 18.21736196 | 82.41637409 | 60.38082896 | 38.01325118 | 54.23805064 | 32.51346997 | 65.72597036 | 47.39000000 | 43.47626753 | ... | 28.34650434 | 40.50305117 | 77.30967665 | 41.40467767 | 23.85318717 | 23.68566393 | 48.15052124 | 79.56803513 | 27.88706274 | 30.01295264 |
| 2013-08-07 | 31.84516887 | 18.16921883 | 81.53454465 | 60.34578796 | 37.75068193 | 54.31202002 | 32.36035945 | 65.50164964 | 47.10000000 | 43.23874037 | ... | 28.19288924 | 40.52972131 | 77.19980174 | 41.22964615 | 23.61393755 | 23.19201856 | 48.07243931 | 79.17498438 | 27.57383161 | 30.11900548 |
| 2013-08-08 | 31.54928679 | 18.27513373 | 80.90042011 | 60.22638901 | 38.16571074 | 55.11643707 | 32.35135295 | 65.48370398 | 47.51000000 | 43.19475386 | ... | 28.02120177 | 40.52083127 | 77.57168605 | 41.57970919 | 23.87711213 | 23.26653107 | 48.21558951 | 79.84604489 | 28.01994868 | 30.11900548 |
| 2013-08-09 | 31.80388300 | 17.90924591 | 82.40646589 | 59.36939001 | 37.86926159 | 55.00548300 | 32.32433344 | 65.97720956 | 47.18000000 | 43.10678084 | ... | 27.86758667 | 40.27190997 | 77.20825366 | 41.39495370 | 23.99673694 | 23.26653107 | 48.41079433 | 79.15581519 | 27.99147312 | 29.80084697 |
| 2013-08-12 | 31.96214550 | 18.12107570 | 81.69307578 | 61.05595360 | 38.14877079 | 54.24729681 | 32.33333995 | 65.31322023 | 47.20000000 | 43.46747023 | ... | 27.76818866 | 40.35192038 | 76.50187303 | 41.53108932 | 24.28383649 | 23.12682011 | 48.45634212 | 78.90656401 | 28.10537535 | 29.24165929 |
| ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
| 2017-05-19 | 55.50327007 | 44.83282860 | 151.06072036 | 150.70113045 | 63.42995100 | 87.59994036 | 42.31486674 | 118.75601310 | 136.43000000 | 79.14250576 | ... | 40.43133035 | 65.31985198 | 78.78898294 | 61.13598414 | 26.81497802 | 51.24320268 | 68.90686651 | 116.46112494 | 39.54646594 | 59.92967369 |
| 2017-05-22 | 55.45382835 | 45.81435227 | 146.97179877 | 151.61679784 | 63.29454092 | 87.91434122 | 42.87370534 | 120.45421034 | 138.86000000 | 79.97056004 | ... | 40.86051697 | 66.15263846 | 79.13518133 | 62.09836493 | 26.73836380 | 51.49936943 | 69.83126290 | 116.57989212 | 39.65494752 | 59.92967369 |
| 2017-05-23 | 58.00502088 | 46.26049939 | 140.27992953 | 151.42972601 | 63.68142686 | 87.62941544 | 42.82468441 | 119.90450488 | 139.52000000 | 79.84391644 | ... | 41.31896631 | 62.67453024 | 79.41406336 | 62.65396622 | 26.62344247 | 52.19890171 | 69.74275686 | 116.59968666 | 40.76934918 | 61.04261076 |
| 2017-05-24 | 58.56865644 | 46.36955757 | 132.66057320 | 150.97681525 | 63.76847620 | 88.39576754 | 42.67762162 | 119.70818149 | 141.12000000 | 79.99978549 | ... | 41.61159355 | 63.13501218 | 79.13518133 | 62.23726525 | 26.89159225 | 52.01106475 | 70.75565928 | 117.87643393 | 40.13818363 | 61.90712437 |
| 2017-05-25 | 58.63787484 | 47.60885514 | 131.61341035 | 151.49864721 | 64.14569000 | 89.51582062 | 43.08939743 | 120.83704093 | 142.85000000 | 80.21410542 | ... | 42.29439045 | 64.10496348 | 78.61588375 | 62.57459460 | 26.77667091 | 51.53652928 | 70.93267135 | 118.07437924 | 40.31569894 | 62.18535864 |
| 2017-05-26 | 58.84553005 | 48.32269053 | 133.79749286 | 151.24265417 | 63.89421413 | 89.40774532 | 43.83451556 | 120.63090138 | 141.89000000 | 80.67197073 | ... | 42.23586500 | 64.51645797 | 78.42355131 | 62.22734380 | 26.81497802 | 50.82472608 | 70.89333534 | 118.00509838 | 39.89163460 | 62.21516946 |
| 2017-05-30 | 59.69592756 | 47.54936885 | 132.63065426 | 151.30172949 | 63.85552554 | 89.43722040 | 44.11883696 | 121.49472426 | 142.41000000 | 82.61059193 | ... | 42.33340741 | 64.38909063 | 77.99080333 | 62.12812929 | 27.12143491 | 51.20039999 | 71.20802347 | 117.21331713 | 39.31964082 | 61.86737662 |
| 2017-05-31 | 59.66626253 | 47.99551598 | 133.26892494 | 150.40575387 | 63.85552554 | 90.16427240 | 44.76591323 | 122.18185609 | 141.86000000 | 83.54580618 | ... | 42.61628041 | 65.35904193 | 77.41380602 | 63.02105992 | 27.08312780 | 51.54641544 | 71.43420556 | 117.98530385 | 39.51688006 | 61.88725050 |
| 2017-06-01 | 60.05190791 | 48.63003633 | 136.72954883 | 150.81928108 | 64.52290380 | 91.51030109 | 45.19729742 | 122.98678196 | 141.38000000 | 80.08746182 | ... | 42.54800072 | 65.30025701 | 77.60613845 | 63.55681830 | 27.19804914 | 51.84300011 | 72.60445204 | 121.40975776 | 39.99025421 | 62.24498027 |
| 2017-06-02 | 60.13101466 | 49.09601221 | 137.42765739 | 153.05429719 | 65.04519983 | 91.94260228 | 45.58946486 | 123.42850956 | 143.48000000 | 78.82102586 | ... | 42.21635651 | 65.52559923 | 76.45214383 | 63.94375491 | 27.12143491 | 52.39662482 | 72.76179611 | 122.66671050 | 39.54646594 | 62.10586314 |
| 2017-06-05 | 59.72559259 | 49.31412858 | 135.20368297 | 151.55772252 | 65.29667569 | 91.68715157 | 45.70711509 | 124.25306776 | 143.59000000 | 76.71679380 | ... | 41.72864445 | 65.79013140 | 77.04837438 | 63.39807508 | 26.73836380 | 52.59434793 | 72.96831019 | 122.65681324 | 39.90149656 | 62.27479108 |
| 2017-06-06 | 59.42894229 | 49.31412858 | 130.94522072 | 152.06970859 | 65.64487304 | 90.08567218 | 45.45220625 | 124.00766354 | 143.03000000 | 78.03193884 | ... | 41.48478841 | 66.24081585 | 78.09658617 | 63.05082428 | 26.81497802 | 52.09015400 | 73.08631824 | 122.89434761 | 39.70425733 | 62.58283617 |
| 2017-06-07 | 59.95302448 | 50.42453920 | 130.26705812 | 152.97553010 | 66.49602213 | 90.32147283 | 45.64828997 | 124.22361925 | 143.62000000 | 79.18147302 | ... | 41.45552569 | 66.49555053 | 77.80808751 | 63.10043153 | 26.96820647 | 52.85138798 | 73.02731422 | 123.07249839 | 39.90149656 | 62.86107043 |
| 2017-06-08 | 59.47838401 | 50.98965889 | 125.57975776 | 152.60138644 | 66.50569428 | 89.96777186 | 45.80515695 | 123.85060483 | 142.63000000 | 80.75863709 | ... | 41.18240693 | 66.69150029 | 77.52920548 | 62.95160976 | 26.81497802 | 52.97002185 | 72.74212810 | 123.88407418 | 40.81865898 | 62.18535864 |
| 2017-06-09 | 58.54887976 | 49.83959075 | 128.01316475 | 146.68400898 | 67.38585980 | 90.48849829 | 46.36399555 | 123.50703891 | 138.05000000 | 76.99695384 | ... | 41.52380538 | 64.07557102 | 78.98131538 | 62.84247379 | 26.58513535 | 53.35558192 | 71.88656975 | 123.72571793 | 41.75554533 | 62.19529558 |
| 2017-06-12 | 58.33133621 | 49.05635469 | 130.59616644 | 143.17887358 | 67.25044972 | 90.74394899 | 46.24634532 | 123.97821503 | 137.25000000 | 78.11370356 | ... | 41.13363573 | 62.92926493 | 79.75064513 | 62.25710816 | 27.04482069 | 53.40501269 | 70.71632326 | 123.71582066 | 42.33740107 | 61.45996216 |
| 2017-06-13 | 58.61809816 | 49.02661155 | 131.27432905 | 144.33084223 | 67.38585980 | 91.37275071 | 46.53066671 | 124.73406004 | 139.09000000 | 79.57331502 | ... | 41.54331386 | 63.46812677 | 79.77949499 | 62.78294509 | 26.85328513 | 53.17763111 | 71.46370757 | 124.18099215 | 42.53464030 | 61.67360633 |
| 2017-06-14 | 58.71698159 | 48.96712526 | 130.23713918 | 142.92288054 | 68.20799244 | 92.25700314 | 46.70714206 | 124.91075109 | 138.25000000 | 79.28922957 | ... | 41.97472897 | 63.56610164 | 78.92361565 | 63.22941040 | 26.54682824 | 53.04911109 | 71.82756572 | 124.30965660 | 42.63325991 | 61.94235833 |
| 2017-06-15 | 58.54887976 | 48.68952261 | 130.79562603 | 142.06628846 | 68.28536963 | 92.77772957 | 47.17774299 | 124.69479537 | 137.52000000 | 78.12349961 | ... | 42.63165652 | 63.54650667 | 79.10633146 | 62.81270944 | 26.61386569 | 53.34569576 | 71.40470355 | 124.54719098 | 42.27822930 | 62.10161877 |
| 2017-06-16 | 58.84553005 | 48.37226243 | 129.80830106 | 140.07741950 | 68.72061632 | 90.91097444 | 47.26598066 | 125.21505233 | 137.84000000 | 78.40758506 | ... | 43.18073029 | 63.42893681 | 80.28917595 | 63.19964605 | 27.29381692 | 53.43467116 | 71.57188162 | 124.62636910 | 42.49519245 | 62.26087921 |
| 2017-06-19 | 59.87391774 | 49.23481354 | 129.24981421 | 144.08469508 | 69.00110863 | 92.10962773 | 47.93266531 | 125.42119188 | 140.35000000 | 78.73085471 | ... | 43.23955963 | 64.56544541 | 79.58716256 | 63.47744669 | 27.57154347 | 53.55330503 | 72.70279208 | 125.78434918 | 42.76146542 | 62.73866054 |
| 2017-06-20 | 59.65637419 | 47.61876952 | 123.24608056 | 142.77519225 | 68.88504285 | 91.61837639 | 47.81501508 | 124.20398692 | 140.91000000 | 77.58471685 | ... | 43.33760851 | 63.93840619 | 79.15441457 | 62.85239525 | 27.13101169 | 53.26660652 | 72.68312408 | 125.91301364 | 42.36698695 | 62.70879921 |
| 2017-06-21 | 59.12240366 | 48.01534474 | 119.84529455 | 143.62193844 | 69.00110863 | 93.94690777 | 47.61893136 | 124.77332472 | 144.24000000 | 78.34880876 | ... | 43.15131563 | 64.78099015 | 78.31776847 | 63.26909621 | 26.70005669 | 52.93047722 | 73.16499028 | 127.43719255 | 41.74568337 | 62.70879921 |
| 2017-06-22 | 59.93324780 | 48.55072128 | 120.43399427 | 143.38563718 | 70.78078399 | 94.68378480 | 48.30522438 | 119.83579169 | 143.69000000 | 79.66147947 | ... | 42.42575385 | 65.23167459 | 77.97157008 | 63.32862492 | 26.78624769 | 53.23694805 | 73.30266633 | 128.29986262 | 41.71609749 | 63.21644187 |
| 2017-06-23 | 59.10262697 | 48.21363235 | 119.47610998 | 144.02561977 | 70.25848796 | 94.14340831 | 48.11894484 | 120.48365885 | 145.41000000 | 79.88678863 | ... | 42.56302230 | 66.16243594 | 78.48125104 | 63.33854637 | 27.23635625 | 53.71148352 | 73.57801845 | 128.00239018 | 41.35120491 | 62.48981610 |
| 2017-06-26 | 58.57854478 | 48.36234805 | 121.52159207 | 143.57270901 | 70.35520945 | 94.31043377 | 47.95227368 | 120.09101209 | 144.96000000 | 78.92677572 | ... | 42.76892496 | 65.99587865 | 78.12543603 | 63.56673975 | 27.95461459 | 54.05749897 | 73.49934641 | 127.97264293 | 41.75554533 | 62.43009343 |
| 2017-06-27 | 58.22256443 | 48.08474540 | 121.69121741 | 141.51491885 | 70.01668424 | 93.85848253 | 47.71697322 | 119.94376955 | 142.54000000 | 76.54633554 | ... | 43.14151074 | 63.78164638 | 78.00041995 | 63.92391201 | 27.75350225 | 53.87954816 | 72.74212810 | 127.16946735 | 41.95278457 | 62.46990854 |
| 2017-06-28 | 58.73675827 | 48.82832394 | 116.45278767 | 143.58255490 | 70.52930812 | 94.69360982 | 47.53069368 | 121.46527575 | 143.81000000 | 77.58471685 | ... | 43.30819385 | 64.67321778 | 78.40431807 | 64.82428373 | 28.28980181 | 54.34419748 | 72.91914017 | 127.42727680 | 42.37684891 | 62.65903032 |
| 2017-06-29 | 58.27398382 | 49.19515602 | 115.79424221 | 141.46568942 | 70.10373358 | 94.08445815 | 47.77579833 | 120.72906307 | 141.24000000 | 76.15449354 | ... | 43.27877918 | 62.88027749 | 77.60613845 | 64.10898129 | 28.12560699 | 54.27499439 | 72.23075989 | 126.81250043 | 43.38276899 | 62.21111032 |
| 2017-06-30 | 58.77942143 | 49.88916265 | 116.33305213 | 141.80044954 | 70.13275003 | 92.87597984 | 47.65814810 | 121.40637874 | 141.44000000 | 76.21326984 | ... | 42.94541296 | 63.01744232 | 77.63498832 | 64.41695873 | 27.74892476 | 54.79896064 | 72.53561401 | 127.31820357 | 43.30387330 | 62.09166499 |
1009 rows × 519 columns
Stock Example¶
Let's see what a single stock looks like from the closing prices. For this example and future display examples in this project, we'll use Apple's stock (AAPL). If we tried to graph all the stocks, it would be too much information.
In [5]:
apple_ticker = 'AAPL'
project_helper.plot_stock(close[apple_ticker], '{} Stock'.format(apple_ticker))
The Alpha Research Process¶
In this project you will code and evaluate a "breakout" signal. It is important to understand where these steps fit in the alpha research workflow. The signal-to-noise ratio in trading signals is very low and, as such, it is very easy to fall into the trap of overfitting to noise. It is therefore inadvisable to jump right into signal coding. To help mitigate overfitting, it is best to start with a general observation and hypothesis; i.e., you should be able to answer the following question before you touch any data:
What feature of markets or investor behaviour would lead to a persistent anomaly that my signal will try to use?
Ideally the assumptions behind the hypothesis will be testable before you actually code and evaluate the signal itself. The workflow therefore is as follows:
In this project, we assume that the first three steps area done ("observe & research", "form hypothesis", "validate hypothesis"). The hypothesis you'll be using for this project is the following:
- In the absence of news or significant investor trading interest, stocks oscillate in a range.
- Traders seek to capitalize on this range-bound behaviour periodically by selling/shorting at the top of the range and buying/covering at the bottom of the range. This behaviour reinforces the existence of the range.
- When stocks break out of the range, due to, e.g., a significant news release or from market pressure from a large investor:
- the liquidity traders who have been providing liquidity at the bounds of the range seek to cover their positions to mitigate losses, thus magnifying the move out of the range, and
- the move out of the range attracts other investor interest; these investors, due to the behavioural bias of herding (e.g., Herd Behavior) build positions which favor continuation of the trend.
Using this hypothesis, let start coding..
Compute the Highs and Lows in a Window¶
You'll use the price highs and lows as an indicator for the breakout strategy. In this section, implement get_high_lows_lookback to get the maximum high price and minimum low price over a window of days. The variable lookback_days contains the number of days to look in the past. Make sure this doesn't include the current day.
In [6]:
def get_high_lows_lookback(high, low, lookback_days):
"""
Get the highs and lows in a lookback window.
Parameters
----------
high : DataFrame
High price for each ticker and date
low : DataFrame
Low price for each ticker and date
lookback_days : int
The number of days to look back
Returns
-------
lookback_high : DataFrame
Lookback high price for each ticker and date
lookback_low : DataFrame
Lookback low price for each ticker and date
"""
#TODO: Implement function
# See Volatility/ rolling Windows
# Calculate the rolling maximum for the given lookback_days, then shift to exclude the current day
lookback_high = high.rolling(window=lookback_days).max().shift(1)
# Calculate the rolling minimum for the given lookback_days, then shift to exclude the current day
lookback_low = low.rolling(window = lookback_days).min().shift(1)
return lookback_high, lookback_low
project_tests.test_get_high_lows_lookback(get_high_lows_lookback)
Tests Passed
View Data¶
Let's use your implementation of get_high_lows_lookback to get the highs and lows for the past 50 days and compare it to it their respective stock. Just like last time, we'll use Apple's stock as the example to look at.
In [7]:
lookback_days = 50
lookback_high, lookback_low = get_high_lows_lookback(high, low, lookback_days)
project_helper.plot_high_low(
close[apple_ticker],
lookback_high[apple_ticker],
lookback_low[apple_ticker],
'High and Low of {} Stock'.format(apple_ticker))
Compute Long and Short Signals¶
Using the generated indicator of highs and lows, create long and short signals using a breakout strategy. Implement get_long_short to generate the following signals:
| Signal | Condition |
|---|---|
| -1 | Low > Close Price |
| 1 | High < Close Price |
| 0 | Otherwise |
In this chart, Close Price is the close parameter. Low and High are the values generated from get_high_lows_lookback, the lookback_high and lookback_low parameters.
In [8]:
def get_long_short(close, lookback_high, lookback_low):
"""
Generate the signals long, short, and do nothing.
Parameters
----------
close : DataFrame
Close price for each ticker and date
lookback_high : DataFrame
Lookback high price for each ticker and date
lookback_low : DataFrame
Lookback low price for each ticker and date
Returns
-------
long_short : DataFrame
The long, short, and do nothing signals for each ticker and date
"""
# Define the conditions for the long and short signals
conditions = [
lookback_low > close, # condition for -1
lookback_high < close, # condition for 1
(lookback_high >= close) & (lookback_low <= close) # condition for 0
]
# Define the corresponding values for the conditions
choices = [-1, 1, 0]
# Use numpy's select method to apply the conditions and get the signals
long_short = np.select(conditions, choices, default=0)
# Convert the numpy array to a DataFrame
long_short = pd.DataFrame(long_short, index=close.index, columns=close.columns)
return long_short
project_tests.test_get_long_short(get_long_short)
Tests Passed
View Data¶
Let's compare the signals you generated against the close prices. This chart will show a lot of signals. Too many in fact. We'll talk about filtering the redundant signals in the next problem.
In [9]:
signal = get_long_short(close, lookback_high, lookback_low)
project_helper.plot_signal(
close[apple_ticker],
signal[apple_ticker],
'Long and Short of {} Stock'.format(apple_ticker))
Filter Signals¶
That was a lot of repeated signals! If we're already shorting a stock, having an additional signal to short a stock isn't helpful for this strategy. This also applies to additional long signals when the last signal was long.
Implement filter_signals to filter out repeated long or short signals within the lookahead_days. If the previous signal was the same, change the signal to 0 (do nothing signal). For example, say you have a single stock time series that is
[1, 0, 1, 0, 1, 0, -1, -1]
Running filter_signals with a lookahead of 3 days should turn those signals into
[1, 0, 0, 0, 1, 0, -1, 0]
To help you implement the function, we have provided you with the clear_signals function. This will remove all signals within a window after the last signal. For example, say you're using a windows size of 3 with clear_signals. It would turn the Series of long signals
[0, 1, 0, 0, 1, 1, 0, 1, 0]
into
[0, 1, 0, 0, 0, 1, 0, 0, 0]
clear_signals only takes a Series of the same type of signals, where 1 is the signal and 0 is no signal. It can't take a mix of long and short signals. Using this function, implement filter_signals.
For implementing filter_signals, we don't reccommend you try to find a vectorized solution. Instead, you should use the iterrows over each column.
In [10]:
def clear_signals(signals, window_size):
"""
Clear out signals in a Series of just long or short signals.
Remove the number of signals down to 1 within the window size time period.
Parameters
----------
signals : Pandas Series
The long, short, or do nothing signals
window_size : int
The number of days to have a single signal
Returns
-------
signals : Pandas Series
Signals with the signals removed from the window size
"""
# Start with buffer of window size
# This handles the edge case of calculating past_signal in the beginning
clean_signals = [0]*window_size
for signal_i, current_signal in enumerate(signals):
# Check if there was a signal in the past window_size of days
has_past_signal = bool(sum(clean_signals[signal_i:signal_i+window_size]))
# Use the current signal if there's no past signal, else 0/False
clean_signals.append(not has_past_signal and current_signal)
# Remove buffer
clean_signals = clean_signals[window_size:]
# Return the signals as a Series of Ints
return pd.Series(np.array(clean_signals).astype(np.int), signals.index)
def filter_signals(signal, lookahead_days):
"""
Filter out signals in a DataFrame.
Parameters
----------
signal : DataFrame
The long, short, and do nothing signals for each ticker and date
lookahead_days : int
The number of days to look ahead
Returns
-------
filtered_signal : DataFrame
The filtered long, short, and do nothing signals for each ticker and date
"""
#TODO: Implement function
# Create an empty DataFrame to store the filtered signals
filtered_signal = pd.DataFrame(index=signal.index)
# Iterate over each column (ticker)
for ticker, col_data in signal.iteritems():
# Separate out the long and short signals
long_signals = col_data == 1
short_signals = col_data == -1
# Use clear_signals function to filter out repeated signals
cleaned_long_signals = clear_signals(long_signals, lookahead_days)
cleaned_short_signals = clear_signals(short_signals, lookahead_days)
# Combine the cleaned long and short signals
filtered_signal[ticker] = cleaned_long_signals - cleaned_short_signals
# Return as int64 type
return filtered_signal.astype(np.int64)
project_tests.test_filter_signals(filter_signals)
Tests Passed
View Data¶
Let's view the same chart as before, but with the redundant signals removed.
In [11]:
signal_5 = filter_signals(signal, 5)
signal_10 = filter_signals(signal, 10)
signal_20 = filter_signals(signal, 20)
for signal_data, signal_days in [(signal_5, 5), (signal_10, 10), (signal_20, 20)]:
project_helper.plot_signal(
close[apple_ticker],
signal_data[apple_ticker],
'Long and Short of {} Stock with {} day signal window'.format(apple_ticker, signal_days))
Lookahead Close Prices¶
With the trading signal done, we can start working on evaluating how many days to short or long the stocks. In this problem, implement get_lookahead_prices to get the close price days ahead in time. You can get the number of days from the variable lookahead_days. We'll use the lookahead prices to calculate future returns in another problem.
In [12]:
def get_lookahead_prices(close, lookahead_days):
"""
Get the lookahead prices for `lookahead_days` number of days.
Parameters
----------
close : DataFrame
Close price for each ticker and date
lookahead_days : int
The number of days to look ahead
Returns
-------
lookahead_prices : DataFrame
The lookahead prices for each ticker and date
"""
# Shift the close prices backwards by `lookahead_days` to get the lookahead prices
lookahead_prices = close.shift(-lookahead_days)
return lookahead_prices
project_tests.test_get_lookahead_prices(get_lookahead_prices)
Tests Passed
View Data¶
Using the get_lookahead_prices function, let's generate lookahead closing prices for 5, 10, and 20 days.
Let's also chart a subsection of a few months of the Apple stock instead of years. This will allow you to view the differences between the 5, 10, and 20 day lookaheads. Otherwise, they will mesh together when looking at a chart that is zoomed out.
In [13]:
lookahead_5 = get_lookahead_prices(close, 5)
lookahead_10 = get_lookahead_prices(close, 10)
lookahead_20 = get_lookahead_prices(close, 20)
project_helper.plot_lookahead_prices(
close[apple_ticker].iloc[150:250],
[
(lookahead_5[apple_ticker].iloc[150:250], 5),
(lookahead_10[apple_ticker].iloc[150:250], 10),
(lookahead_20[apple_ticker].iloc[150:250], 20)],
'5, 10, and 20 day Lookahead Prices for Slice of {} Stock'.format(apple_ticker))
Lookahead Price Returns¶
Implement get_return_lookahead to generate the log price return between the closing price and the lookahead price.
In [14]:
def get_return_lookahead(close, lookahead_prices):
"""
Calculate the log returns from the lookahead days to the signal day.
Parameters
----------
close : DataFrame
Close price for each ticker and date
lookahead_prices : DataFrame
The lookahead prices for each ticker and date
Returns
-------
lookahead_returns : DataFrame
The lookahead log returns for each ticker and date
"""
# Calculate the log returns
lookahead_returns = np.log(lookahead_prices / close)
return lookahead_returns
project_tests.test_get_return_lookahead(get_return_lookahead)
Tests Passed
View Data¶
Using the same lookahead prices and same subsection of the Apple stock from the previous problem, we'll view the lookahead returns.
In order to view price returns on the same chart as the stock, a second y-axis will be added. When viewing this chart, the axis for the price of the stock will be on the left side, like previous charts. The axis for price returns will be located on the right side.
In [15]:
price_return_5 = get_return_lookahead(close, lookahead_5)
price_return_10 = get_return_lookahead(close, lookahead_10)
price_return_20 = get_return_lookahead(close, lookahead_20)
project_helper.plot_price_returns(
close[apple_ticker].iloc[150:250],
[
(price_return_5[apple_ticker].iloc[150:250], 5),
(price_return_10[apple_ticker].iloc[150:250], 10),
(price_return_20[apple_ticker].iloc[150:250], 20)],
'5, 10, and 20 day Lookahead Returns for Slice {} Stock'.format(apple_ticker))
Compute the Signal Return¶
Using the price returns generate the signal returns.
In [16]:
def get_signal_return(signal, lookahead_returns):
"""
Compute the signal returns.
Parameters
----------
signal : DataFrame
The long, short, and do nothing signals for each ticker and date
lookahead_returns : DataFrame
The lookahead log returns for each ticker and date
Returns
-------
signal_return : DataFrame
Signal returns for each ticker and date
"""
# Multiply the signal with the lookahead returns
signal_return = signal * lookahead_returns
return signal_return
project_tests.test_get_signal_return(get_signal_return)
Tests Passed
View Data¶
Let's continue using the previous lookahead prices to view the signal returns. Just like before, the axis for the signal returns is on the right side of the chart.
In [17]:
title_string = '{} day LookaheadSignal Returns for {} Stock'
signal_return_5 = get_signal_return(signal_5, price_return_5)
signal_return_10 = get_signal_return(signal_10, price_return_10)
signal_return_20 = get_signal_return(signal_20, price_return_20)
project_helper.plot_signal_returns(
close[apple_ticker],
[
(signal_return_5[apple_ticker], signal_5[apple_ticker], 5),
(signal_return_10[apple_ticker], signal_10[apple_ticker], 10),
(signal_return_20[apple_ticker], signal_20[apple_ticker], 20)],
[title_string.format(5, apple_ticker), title_string.format(10, apple_ticker), title_string.format(20, apple_ticker)])
In [18]:
project_helper.plot_signal_histograms(
[signal_return_5, signal_return_10, signal_return_20],
'Signal Return',
('5 Days', '10 Days', '20 Days'))
Question: What do the histograms tell you about the signal returns?¶
To answer this question think about the following:
- Are the histograms skewed? If yes, which side?
- Is the skewing because of outliers? What could be the reason for this?
- Which side do the outliers appear and what does it signify?
- What will happen to the distribution of the histograms if the outliers are removed?
Another way to check if a distribution is normal is to create a QQ (quantile-quantile) plot. See this link to learn more about QQ plots.
#TODO: Put Answer In this Cell https://chartio.com/learn/charts/histogram-complete-guide/ .
Are the histograms skewed? If yes, which side?*
- Yes, the histograms are skewed. The increasing lookahead windows (from 5 days to 20 days), the histograms become more negatively skewed. This is due to a larger number of stocks giving larger than expected returns.
Is the skewing because of outliers? What could be the reason for this?*
- Yes, the skewing is because of outliers. The outliers, especially in the 10 and 20 day lookahead windows, are causing the histograms to deviate from a normal distribution. As for the reasons, they could be numerous: current events affecting a particular sector, major company announcements, unexpected earnings reports, etc. These events can lead to sudden and substantial changes in stock prices, creating outlier returns.
Which side do the outliers appear and what does it signify?*
- The outliers appear on the positive side, as evidenced by the "lump" at around +0.25 in the 20-day histogram. This signifies that there are certain stocks giving abnormally high returns. These positive outliers suggest the possibility of certain stocks reacting to specific events or situations that lead to larger than expected positive returns.
What will happen to the distribution of the histograms if the outliers are removed?*
- If the outliers are removed, the histograms, especially for the 10 and 20 day lookahead windows, would resemble a more normal distribution. The distributions would be less skewed, more symmetrical, and the mean might still be around zero but without the major distortions caused by the outliers. This would give a more accurate representation of the typical signal return without the influence of the outlier events.
In summary, the histograms tell you that:
- Signal returns for shorter lookahead windows, like 5 days, are more normally distributed.
- As the lookahead window increases, the presence of outliers (particularly on the positive side) makes the histograms negatively skewed.
- These outliers could be due to specific, non-recurring events affecting certain stocks.
- Without these outliers, the distribution of returns for longer lookahead windows would be more normal, which is crucial to consider when evaluating the effectiveness of a trading signal.
Outliers¶
You might have noticed the outliers in the 10 and 20 day histograms. To better visualize the outliers, let's compare the 5, 10, and 20 day signals returns to normal distributions with the same mean and deviation for each signal return distributions.
In [19]:
project_helper.plot_signal_to_normal_histograms(
[signal_return_5, signal_return_10, signal_return_20],
'Signal Return',
('5 Days', '10 Days', '20 Days'))
Kolmogorov-Smirnov Test¶
While you can see the outliers in the histogram, we need to find the stocks that are causing these outlying returns. We'll use the Kolmogorov-Smirnov Test or KS-Test. This test will be applied to teach ticker's signal returns where a long or short signal exits.
In [20]:
# Filter out returns that don't have a long or short signal.
long_short_signal_returns_5 = signal_return_5[signal_5 != 0].stack()
long_short_signal_returns_10 = signal_return_10[signal_10 != 0].stack()
long_short_signal_returns_20 = signal_return_20[signal_20 != 0].stack()
# Get just ticker and signal return
long_short_signal_returns_5 = long_short_signal_returns_5.reset_index().iloc[:, [1,2]]
long_short_signal_returns_5.columns = ['ticker', 'signal_return']
long_short_signal_returns_10 = long_short_signal_returns_10.reset_index().iloc[:, [1,2]]
long_short_signal_returns_10.columns = ['ticker', 'signal_return']
long_short_signal_returns_20 = long_short_signal_returns_20.reset_index().iloc[:, [1,2]]
long_short_signal_returns_20.columns = ['ticker', 'signal_return']
# View some of the data
long_short_signal_returns_5.head(10)
Out[20]:
| ticker | signal_return | |
|---|---|---|
| 0 | A | 0.00732604 |
| 1 | ABC | 0.01639650 |
| 2 | ADP | 0.00981520 |
| 3 | AGENEN | 0.01447001 |
| 4 | AKAM | 0.04400495 |
| 5 | ALGN | 0.01545561 |
| 6 | ALTAIC | 0.02738546 |
| 7 | APC | 0.00305859 |
| 8 | ARMENA | 0.01549105 |
| 9 | BA | 0.08061297 |
This gives you the data to use in the KS-Test.
Now it's time to implement the function calculate_kstest to use Kolmogorov-Smirnov test (KS test) between a distribution of stock returns (the input dataframe in this case) and each stock's signal returns. Run KS test on a normal distribution against each stock's signal returns. Use scipy.stats.kstest perform the KS test. When calculating the standard deviation of the signal returns, make sure to set the delta degrees of freedom to 0.
For this function, we don't reccommend you try to find a vectorized solution. Instead, you should iterate over the groupby function.
Hint: You should compare the signal return of the individual tickers against a normal distribution whose parameters - mean and standard deviation are computed from all tickers.
In [21]:
from scipy.stats import kstest
def calculate_kstest(long_short_signal_returns):
"""
Calculate the KS-Test against the signal returns with a long or short signal.
Parameters
----------
long_short_signal_returns : DataFrame
The signal returns which have a signal.
This DataFrame contains two columns, "ticker" and "signal_return"
Returns
-------
ks_values : Pandas Series
KS static for all the tickers
p_values : Pandas Series
P value for all the tickers
"""
# Get the mean and standard deviation of the signal returns
mean_return = long_short_signal_returns['signal_return'].mean()
std_return = long_short_signal_returns['signal_return'].std(ddof=0)
# Function to apply kstest for each group
def ks_test_on_group(group):
_, p_value = kstest(group['signal_return'], 'norm', args=(mean_return, std_return))
return pd.Series({'ks_value': _, 'p_value': p_value})
# Apply the ks_test_on_group function to each ticker's data
results = long_short_signal_returns.groupby('ticker').apply(ks_test_on_group)
return results['ks_value'], results['p_value']
project_tests.test_calculate_kstest(calculate_kstest)
Tests Passed
View Data¶
Using the signal returns we created above, let's calculate the ks and p values.
In [22]:
ks_values_5, p_values_5 = calculate_kstest(long_short_signal_returns_5)
ks_values_10, p_values_10 = calculate_kstest(long_short_signal_returns_10)
ks_values_20, p_values_20 = calculate_kstest(long_short_signal_returns_20)
print('ks_values_5')
print(ks_values_5.head(10))
print('p_values_5')
print(p_values_5.head(10))
ks_values_5 ticker A 0.17204716 AAL 0.10741137 AAP 0.19687855 AAPL 0.15576784 ABBV 0.16807629 ABC 0.21422267 ABT 0.21364129 ACN 0.28213325 ADBE 0.24259879 ADI 0.19421352 Name: ks_value, dtype: float64 p_values_5 ticker A 0.18758007 AAL 0.72480842 AAP 0.04516745 AAPL 0.24641150 ABBV 0.24737380 ABC 0.02721895 ABT 0.04844943 ACN 0.00587558 ADBE 0.00915703 ADI 0.09911122 Name: p_value, dtype: float64
Find Outliers¶
With the ks and p values calculate, let's find which symbols are the outliers. Implement the find_outliers function to find the following outliers:
- Symbols that pass the null hypothesis with a p-value less than
pvalue_thresholdAND with a KS value aboveks_threshold.
Note: your function should return symbols that meet both requirements above.
In [23]:
def find_outliers(ks_values, p_values, ks_threshold, pvalue_threshold=0.05):
"""
Find outlying symbols using KS values and P-values
Parameters
----------
ks_values : Pandas Series
KS static for all the tickers
p_values : Pandas Series
P value for all the tickers
ks_threshold : float
The threshold for the KS statistic
pvalue_threshold : float
The threshold for the p-value
Returns
-------
outliers : set of str
Symbols that are outliers
"""
#TODO: Implement function
# Conditions for outliers
condition = (p_values < pvalue_threshold) & (ks_values > ks_threshold)
# Extract the symbols that meet the condition
outlier_symbols = ks_values[condition].index.tolist()
# Convert list to set and return
return set(outlier_symbols)
project_tests.test_find_outliers(find_outliers)
Tests Passed
View Data¶
Using the find_outliers function you implemented, let's see what we found.
In [24]:
ks_threshold = 0.8
outliers_5 = find_outliers(ks_values_5, p_values_5, ks_threshold)
outliers_10 = find_outliers(ks_values_10, p_values_10, ks_threshold)
outliers_20 = find_outliers(ks_values_20, p_values_20, ks_threshold)
outlier_tickers = outliers_5.union(outliers_10).union(outliers_20)
print('{} Outliers Found:\n{}'.format(len(outlier_tickers), ', '.join(list(outlier_tickers))))
24 Outliers Found: LINIFO, ORPHAN, GREIGI, DASYST, TARDA, PULCHE, GESNER, ARMENA, KOLPAK, BAKERI, SYLVES, SPRENG, ALTAIC, PRAEST, VVEDEN, BIFLOR, KAUFMA, HUMILI, TURKES, AGENEN, URUMIE, CLUSIA, SAXATI, SCHREN
Show Significance without Outliers¶
Let's compare the 5, 10, and 20 day signals returns without outliers to normal distributions. Also, let's see how the P-Value has changed with the outliers removed.
In [25]:
good_tickers = list(set(close.columns) - outlier_tickers)
project_helper.plot_signal_to_normal_histograms(
[signal_return_5[good_tickers], signal_return_10[good_tickers], signal_return_20[good_tickers]],
'Signal Return Without Outliers',
('5 Days', '10 Days', '20 Days'))
That's more like it! The returns are closer to a normal distribution. You have finished the research phase of a Breakout Strategy. You can now submit your project.
Submission¶
Now that you're done with the project, it's time to submit it. Click the submit button in the bottom right. One of our reviewers will give you feedback on your project with a pass or not passed grade. You can continue to the next section while you wait for feedback.